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Abstract 

Pairs (Hamiltonian system, Lagrangian distribution) , called dynamical Lagrangian 
distributions, appear naturally in Differential Geometry, Calculus of Variations and 
Rational Mechanics. The basic differential invariants of a dynamical Lagrangian distri- 
bution w.r.t. the action of the group of symplectomorphisms of the ambient symplectic 
manifold are the curvature operator and the curvature form. These invariants can be 
seen as generalizations of the classical curvature tensor in Riemannian Geometry. In 
particular, in terms of these invariants one can localize the focal points along extremals 
of the corresponding variational problems. In the present paper we study the behavior 
of the curvature operator, the curvature form and the focal points of a dynamical La- 
grangian distribution after its reduction by arbitrary first integrals in involution. The 
interesting phenomenon is that the curvature form of so-called monotone increasing 
Lagrangian dynamical distributions, which appear naturally in mechanical systems, 
does not decrease after reduction. It also turns out that the set of focal points to the 
given point w.r.t. the monotone increasing dynamical Lagrangian distribution and the 
corresponding set of focal points w.r.t. its reduction by one integral are alternating 
sets on the corresponding integral curve of the Hamiltonian system of the considered 
dynamical distributions. Moreover, the first focal point corresponding to the reduced 
Lagrangian distribution comes before any focal point related to the original dynamical 
distribution. We illustrate our results on the classical V-body problem. 

Key words: curvature operator and form, focal points, reduction by first integrals, 
curves in Lagrangian Grassmannians. 

1 Introduction 

In the present paper smooth objects are supposed to be C°°. The results remain valid 
for the class with a finite and not large k but we prefer not to specify minimal possible 
k. 
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1.1 Dynamical Lagrangian distributions. Let be a symplectic manifold with 
symplectic form a. Lagrangian distribution T> on. W \s a. smooth vector sub-bundle of 
the tangent bundle TW such that each fiber T)x is a Lagrangian subspace of the linear 
symplectic space TxW , i.e., dvcuDx = ^ dimM^ and crA(^^i;^^2) = for all vi,V2 G T>\. For 
example, as a symplectic manifold one can take the cotangent bundle T*M of a manifold 
M with standard symplectic structure and as a Lagrangian distribution one can take the 
distribution n(M) of tangent spaces to the fibers of T*M, namely, 

n(M), = r,(T;(,)M), (li) 

where vr : T*M ^ M is the canonical projection on the base manifold M. 

Let be a smooth function on W. Denote by Tl the Hamiltonian vector field, corre- 
sponding to the function TC: d7i{-) = a{-,T-L), and by e*''^ the Hamiltonian flow generated 
by Ti. The pair (Tl, D) defines the one-parametric family of Lagrangian distributions 
= (e*^)^!?. The pair (Ji.,T)) will be called dynamical Lagrangian distribution. The 
point Ai = e*^^Ao is called focal to Aq w.r.t. the pair {H,T)) along the integral curve 
t ^ e*^Ao of H, if 

(e*^^\PAonPA, /O. (L2) 

Dynamical Lagrangian distributions appear naturally in Differential Geometry, Calcu- 
lus of Variations and Rational Mechanics. The model example can be described as follows: 



Example 1 On a manifold M for a given smooth function L : TM M, which is convex 
on each fiber, consider the following standard problem of Calculus of Variation with fixed 
endpoints qo and qi and fixed time T: 

^ («(•)) = L{q{t),q{t))dt^ mm (L3) 

q{0)=qo, q{T)=qi. (L4) 
Suppose that the Legendre transform H : T*M i-^ R of the function L, 

= max -L(g,X)), qGM,p£T*M, (L5) 

is well defined and smooth on T*M. We will say that the dynamical Lagrangian distribu- 
tions {H,n{M)) is associated with the problem ((0)) - (fO)) ^ The curve q : [0,r] i-^ M, 
satisfying (|1.4j) . is an extremal of the problem (|1.3j) - p.4|) if and only if there exists an 
integral curve 7 : [0,T] T*M of H such that q{t) = 7r(7(t)) for all < t < T. In this 
case the point 7(r) is focal to 7(0) w.r.t. the pair {H, T>) if and only if qi is conjugate to 
go along the extremal q{-) in the classical variational sense for the problem (|1.3|) - (|1.4|) . □ 

^In the model example the Lagrangian distributions are integrable. For application of one-parametric 
families of non-integrable Lagrangian distributions see 
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The group of symplectomorphisms of W acts naturally on Lagrangian distribution and 
Hamiltonian vector fields, therefore it acts also on dynamical Lagrangian distributions. 
Dynamical Lagrangian distributions have richer geometry w.r.t. this action than just 
Lagrangian distribution. For example, all integrable Lagrangian distributions are locally 
equivalent w.r.t. the action of the group of symplectomorphisms of W, while integrable 
dynamical Lagrangian distributions have functional moduli w.r.t. this action. 

First note that for any two vector fields Y, Z tangent to the distribution D the number 
ax{^[Ti.,Y], depends only on the vectors Y{X), Z{X)? Therefore for a given dynamical 

Lagrangian distribution (Ti.,T)) the following bilinear form Q^^'^\-, ■) is defined on each 

yv,w€Vx: Q^^''^\v,w) =ax{[n,Y],Z), Y{X) =v,Z{X) = w. (L6) 

Moreover, from the fact that all Vx are Lagrangian it follows that the form Q^^'^^ is 
symmetric. 

A dynamical Lagrangian distribution (TC, T>) is called regular, if the quadratic forms 

are non-degenerated for any A. A dynamical Lagrangian distribu- 
tions is called monotone (non-decreasing or non-increasing), if the quadratic forms v i— > 

Q\ ' {v,v) are sign-definite (non-negative or non-positive definite) for any A . The regu- 
lar dynamical Lagrangian distributions is called monotone increasing (decreasing), if the 

quadratic forms v ^ Q^^'^\v,v) are positive (negative) definite for any A . 

Remark 1 If 1^ = T*M and V = U{M) are as in (|rT|) . then the form v ^ Q^'^^^'^^ {v, v) 
coincides with the second differential at A of the restriction 7Y|^, m Hamiltonian 

ir(A) 

7i to the fiber T*^y^M. Therefore in this case the dynamical Lagrangian distribution 
(Ti, n(M)) is monotone increasing if and only if the restrictions of Ti on each fiber of T*M 
are strongly convex. Consequently the dynamical Lagrangian distributions (H, n(M)) as- 
sociated with the problem Hl.c{() - H1.4|) is monotone increasing if and only if the restrictions 
of the function L : TM i-^ R on each fiber of TM are strongly convex. □ 

It turns out that under some non-restrictive assumptions on the dynamical Lagrangian 
distribution [Ti, D) (in particular, if this dynamical Lagrangian distribution is regular) 

one can assign to it a special linear operator R^^'^^ on each linear spaces Vx- This 
operator is called the curvature operator of (Ti., D) at A and it is the basic differential 
invariant of dynamical Lagrangian distribution (7Y, V) w.r.t. the action of the group of 
symplectomorphisms of W. Moreover, the following bilinear form 

rf^^\v,w) = Qf^'^\Rf''^\,w), v,weVx (1.7) 

is symmetric. The corresponding quadratic form is called the curvature form of the pair 
{1-L,T>). Besides, the trace of the curvature operator 

pf» = tri?f>) (1.8) 



Here [wi, W2] is the Lie bracket of the vector fields wi and W2, [wi, ^2] = fi o i;2 — i;2 o ui . 
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is called the generalized Ricci curvature of {7i,'D) at A. All these invariants where intro- 
duced in (see also j3] and section 2 below) and the effective method for their calculations 
is given in the recent work j^. Below we present the results of these calculations on sev- 
eral important examples. In all these examples W = T*M for some manifold M and 
P = n(M), a smooth function L : TM i-^ M is given, the functional is as in ^l.'A^ . 

and H : T*M M is as in ([T3]) . 

n 

Example 2 (Natural mechanical system) M = ii", = M" x M", cj = J2 dpi /\ dqi, 

i=l 

= (IK",0), L{q,X) = ^\\X\\^ - U{q) (in this case the function A{q{-)) is the Action 
functional of the natural mechanical system with potential energy U{q)). Then 

Vl<^,,<n: rlffid,.0,^) = ^^i^)- (1-9) 

In other words, in this case the curvature operator can be identified with the Hessian of 
the potential U. 

Example 3 (Riemannian manifold) Let a Riemannian metric G is given on a manifold 
M by choosing an inner product Gq{-,-) on each subspaces TqM for any q ^ M smoothly 
w.r.t. q. Let L{q,X) = ^Gq{X, X). The inner product Gq{-,-) defines the canonical 
isomorphism between T*M and TqM. For any q £ M and p € T*M we will denote by p^ 
the image of p under this isomorphism, namely, the vector € TqM, satisfying 

p{.)=Gq{p\-) (1.10) 

(the operation ^ corresponds to the operation of raising of indexes in the corresponding 
coordinates of co- vectors and vectors). Since the fibers of T*M are linear spaces, one can 
identify Vx (= TxT*^^^M) with T*^^^M, i.e., the operation 1^ is defined also on each Vx 
with values in T^^^x)^- It turns out (see [2]) that 

yv€Vx: {R['^-''^vy = R^{X\v^)X^ (1.11) 

where R^ is the Riemannian curvature tensor of the metric G. The right-hand side of 
(|l.llj) appears in the classical Jacobi equation for Jacobi vector fields along the Riemannian 

geodesies. Also, :;^irR\ ' is exactly the Ricci curvature calculated at A' . Besides, using 
(jl.llj) , the Riemannian curvature tensor R^ can be recovered uniquely from the curvature 

operator R\^ ' . Therefore studying differential invariants of the appropriate integrable 
dynamical Lagrangian distributions, one can obtain the classical Riemannian tensor. 

Example 4 (Mechanical system on a Riemannian manifold) Let G be the metric of the 
previous example and L{q, X) = ^Gq{X, X) — U{q) (in this case the function ^(g(-)) is the 
Action functional of the mechanical system on the Riemannian manifold with potential 
U{q)). Let the operation ^ be as in (|l.lfl|) . Then the curvature operator satisfies 

^veVx: {Rf'^'^'vy = R^{X\v^)X^ + V,T (grade U)(7r(A)), (1.12) 
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where grad^t/ is the gradient of the function U w.r.t. the metric G, i.e., grad^t/ = dU\ 
and V is the Riemannian covariant derivative. □ 

Remark 2 According to Remark ^ the dynamical Lagrangian distributions from Exam- 
ples 2-4 are monotone increasing. □ 

The generalization of different kinds of Riemannian curvatures, using the notion of the 
curvature operator of dynamical Lagrangian distributions, leads to the generalization of 
several classical results of Riemannian geometry. In 7 for the given monotone increasing 
dynamical Lagrangian distribution (Ti.,!)) the estimates of intervals between two consec- 
utive focal points w.r.t. the pair {71, T)) along the integral curve "y{t) = e*^Ao of 7i were 
obtained in terms of the curvature form of the pair (TtjV). This result is the general- 
ization of the classical Ranch Comparison Theorem in Riemannian geometry, which gives 
the lower and upper bounds of the interval between consecutive conjugate points along 
the Riemannian geodesies in terms of upper bound for the sectional curvatures and lower 
bound for the Ricci curvature respectively. In recent work fT it was shown that the Hamil- 
tonian flow, generated by a vector field Ti. on the compact level set of 7Y, is hyperbolic, if 
there exists a Lagrangian distribution D such that the dynamical Lagrangian distribution 
(Ti.j'D) is monotone (increasing or decreasing) and the curvature form of so-called reduc- 
tion of this dynamical distribution by Hamiltonian TC on this level set is negative definite. 
This is an analog of the classical theorem about hyperbolicity of geodesic flows of negative 
sectional curvature on a compact Riemannian manifold. 

1.2 The reduction by the first integrals. The subject of the present paper is the 
behavior of the curvature form and the focal points after the reduction of the dynamical 
Lagrangian distribution {7i,T>) by the arbitrary s first integral gi, . . . ,gs in involution of 
the Hamiltonian 7i, i.e., s functions on W such that 

{n,gi} = 0, {gi,gj} = 0, VI < i, j < s (1.13) 

(here {h,g} is the Poisson bracket of the functions h and g, {h,g} = dg{h)). This problem 
appears naturally in the framework of mechanical systems and variational problems with 
symmetries. Let Q = {gi, . . . , gs) and 

s 

Vl = (p|ker dA5i) nPA+span(5i(A),...,5.(A)). (1.14) 

Obviously, is a Lagrangian distribution. The pair {Ti, D^) is called the reduction by the 
tuple G of s first integrals of Ti in involution or shortly the Q-reduction of the dynamical 
Lagrangian distribution (Ti,!)). The following example justifies the word "reduction" in 
the previous definition: 

Example 5 Assume that we have one first integral g of Ti. such that the Hamiltonian 
vector field g, corresponding to the first integral g, preserves the distribution Dx, namely, 

ie'^),V = V. (1.15) 
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Fixing some value c of g, one can define (at least locally) the following quotient manifold: 



Wg,, = g~\c)/C, 

where C is the line foliation of the integral curves of the vector field g. The symplectic 
form aoiW induces the symplectic form on a manifold Wg^c-, making it symplectic too. 
Besides, if we denote by <1> : g~^{c) ^ Wg^c the canonical projection on the quotient set, 
the vector field $*(7Y) is well defined Hamiltonian vector field on Wg^c, because by our 
assumptions the vector fields Ti. and g commute. Actually we have described the standard 
reduction of the Hamiltonian systems on the level set of the first integral, commonly 
used in Mechanics. In addition, by (|1.15|) . <I>*(I'^) is well defined Lagrangian distribution 
on Wg^c- So, to any dynamical Lagrangian distribution {'H,'D) on W one can associate 
the dynamical Lagrangian distribution {^i^TC, ^^I?^) on the symplectic manifold Wg^c of 
smaller dimension. It turns out (see subsection 2.2 below) that the curvature form of the 
g-reduction {H,V^) at A € g^^{c) is equal to the pull-back by $ of the curvature form 
of the dynamical Lagrangian distribution (<I>*7Y, $^,1?^). So, instead of (7Y, P^) one can 
work with (<I>*'H, ^^V^) on the reduced symplectic space Wg^c- This is the essence of the 
reduction on the level set of the first integral. □ 

Remark 3 Suppose now that W = T*M for some manifold M and V = n(M). In this 
case if is a first integral of TC, which is "linear w.r.t. the impulses", i.e., there exists a 
vector field V on M such that 

g{p,q)=piV{q)), qeM,peT;M, (1.16) 

then it satisfies ()1.15|) . If we denote by V the line foliation of integral curves of V, then 
the reduced symplectic space Wg^c can be identified with T*(M/V) and the distribution 
n(M)^ can be identified with n(M/V). So, after reduction we work with the dynamical 
Lagrangian distribution (<5*7Y, n(M/V)) on the reduced phase space T*{M/V) instead of 
(7?,n(M)). □ 

In view of the previous example the following analog of the notion of the focal points 
along the extremal w.r.t. the ^-reduction of the pair (7Y, V) is natural: The point Ai = 
e*^^Ao is called focal to Aq w.r.t. the ^-reduction of the pair {Ti,!)) along the integral 
curve t I— > e*^Ao oiTi, if 

((e*^«)Xn2?l)/span(5i(A),...,5s(A)) /O. (1.17) 

In the situation, described in Example El the point Ai = e*^^Ao is focal to Aq w.r.t. the 
^'-reduction of the pair [H^V) along the curve t i— e*^Ao if and only if <I>(Ai) is focal 
to $(Ao) w.r.t. the pair (($)*??, ($)*P^') along the curve t ^ ^>(e*'^Ao) in Wg^c- We 
illustrate the meaning of the focal points of the reduction from the variational point of 
view on the following two examples. In both examples W = T*M for some manifold 
M and T) = n(M), a fiber-wise convex and smooth function L : TM i-^ M is given and 
H : T*M M is as in ((ESI). 
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Example 6 Assume that the Hamiltonian H admits a first integral g, satisfying H1.16() . 
It is well known that g, satisfying ()1.16() . is the first integral of H if and only if the flow e*^ 
induces the one-parametric family of fiber- wise diffeomorphisms on TM, which preserve 
the function L, i.e., L o (e*^)* = L. 

Let Vi : M 1-^ M be an integral curve of V and a(-) be a function on Vi such that 

a{Vi{s))=s, sGM. 

Fix some real c. Then for the given point qo, and the time T consider the following 
variational problem 

[ L{q{t),q{t))dt-ca{q{T)) ^ mm, (1.18) 
Jo 

g(0) = qo, q{T) = Vi. (1.19) 

The curve q : [0,T] i— > M, satisfying (|1.19)) . is an extremal of the problem if 
and only if there exists an integral curve 7 : [0, T] 1— > g~^{c) of H, such that q{t) = 7r[-f{t)) 
for all < t < T. In this case the point 7(0) is focal to 7(T) w.r.t. the (7-reduction of the 
pair (H, T>) if and only if the point qo is focal to the point q(T) along the extremal q{-) in 
the classical variational sense for the problem (fTT8|l - (ri~T9] l. □ 

Example 7 Suppose that g = H. For given real c and points qo, qi consider the following 
variational problem with free terminal time 

[ L{q{t),q{t))dt-cT ^ mm, T is free, (1.20) 

g(0) = qo, q(T) = q^. (1.21) 

The curve q : [0,T] 1-^ M, satisfying (fT^ . is an extremal of the problem ((On)) - ((r7T|) if 
and only if there exists an integral curve 7 : [0, T] 1— > H~^{c) of H, such that q{t) = vr(7(t)) 
for all < t < T. In this case the point 7(0) is focal to 7(r) w.r.t. the //-reduction of the 
pair (H,!)) if and only if the point qo is focal to the point qi along the extremal q{-) in 
the classical variational sense for the problem H1.20() - H1.21() . Actually the considered case 
can be seen as a particular case of the previous example. For this one can pass to the 
extended (configuration) space M = M x M instead of M and take the following function 
L : TM ^ R instead of L: 

L(g-X)t/L(<?,-)y, 

y 

where g G M such that q = {q,t), q G M, t G R and X e TgM such that X = {X,y), 
Y = {X,y), X £ TgM, y eTtR = R (it is well known that {t, H) is the pair of conjugate 
variables for function L, so as the field V one takes ^). □ 

1.3 Description of main results. For the reduced dynamical Lagrangian distribu- 
tion (7i, P^) one can also define the curvature operator R^^'^ ^ and the curvature form 
(W,© ) each linear spaces V^. The natural problem is to find the relation between 
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R^x^'^^ (or r^x^'^^) and their reduced analogs R^^''^ ^ (or r^^''^ ^) on the hnear space 

s 

ker dxQi^ D Vx (which is the intersection of the corresponding spaces of definition T>x 

i=l 

and T)^). We solve this problem in section 2 for regular dynamical distributions. It gives 
an effective and flexible method to compute and evaluate the curvature of Hamiltonian 
systems arising in Rational Mechanics and geometric variational problems. The interesting 
phenomenon is that the curvature form of a monotone increasing Lagrangian dynamical 
distribution does not decrease after reduction. More precisely, for such distribution the 
quadratic form 



^'^''''\v,v)-rf'''\v,v), V G (f|ker dxg^) nVx 



i=l 



is always non-negative definite of rank not greater than s, where s is the number of the 
first integrals in the tuple Q. 

Further, in section 3 we show that the set of focal points to the given point along 
an integral curve w.r.t. the monotone increasing (or decreasing) dynamical Lagrangian 
distribution and the corresponding set of focal points w.r.t. its reduction by one integral 
are alternating sets on the curve and the first focal point w.r.t. the reduction comes before 
any focal point w.r.t. the original dynamical Lagrangian distribution. In view of Examples 
ini and [7| this result looks natural: The reduction enlarge the set of admissible curves in 
the corresponding variational problems (instead of the problem with fixed endpoints ant 
terminal time one obtains the problem with variable endpoints or free terminal time). 
This justifies the fact that the first focal point of the reduction comes sooner. Besides, 
for the mentioned examples the last fact and the alternation of focal points are also a 
consequence of the classical Courant Minimax Principle, applied to the second variation 
along the reference extremal in the corresponding variational problems. 

In addition, we demonstrate our results on the classical N - body problem. 



2 Curvature and reduction 



2.1 Curvature operator and curvature form. For the construction of the curva- 
ture operator of the dynamical Lagrangian distribution we use the theory of curves in the 
Lagrange Grassmannian, developed in [2] and 0. The curve 

t^Jx{t)'^^e:'^{V^,^^). (2.1) 

is called the Jacobi curve of the curve t i— > e*^A attached at the point A (w.r.t. the 
dynamical distribution (Ti,!))). It is the curve in the Lagrange Grassmannian L(TxW) 
of the linear symplectic space TxW . Actually, the Jacobi curve is a generalization of the 
space of "Jacobi fields" along the extremal of variational problem of type H1.3() - H1.4() . Note 
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that if A = e^'^X then by (HHJ) we have 

J-x{t) = ef^Ut-t). 

In other words, the Jacobi curves of the same integral curve of Tl attached at two different 
points of this curve are the same, up to symplectic transformation between the correspond- 
ing ambient hnear symplectic spaces and the corresponding shift of the parameterizations. 
Therefore, any differential invariants of the Jacobi curve w.r.t. the action of the linear 
Symplectic group (in other words, any symplectic invariant of the curve) produces the 
invariant of the corresponding dynamical Lagrangian distributions w.r.t. the action of the 
group of symplectomorphisms of the ambient space W. 

Now, following iSI and we describe the construction of the curvature operator of 
the curve in the Lagrange Grassmannian. Let S be 2n-dimensional linear space, endowed 
with symplectic form a. The Lagrange Grassmannian is real analytic manifold. Note 
that the tangent space T\L{T,) to the Lagrangian Grassmannian at the point A can be 
naturally identified with the space of quadratic forms Quad (A) on the linear space A C S. 
Namely, take a curve A{t) £ L(Y,) with A(0) = A. Given some vector / G A, take a curve 

in W such that l{t) G A(t) for all t and /(O) = /. Define the quadratic form 

/-^(|/(0),0. (2.2) 

Using the fact that the spaces A{t) are Lagrangian, it is easy to see that this form depends 
only on J^A(O). So, we have the map from T\L{Y?) to the space Quad(A). A simple 
counting of dimension shows that this mapping is a bijection. 

Remark 4 In the sequel, depending on the context, we will look on the elements of 
T\L{Ti) not only as on the quadratic forms on A(t), but also as on the corresponding 
symmetric bilinear forms on K{t) or on the corresponding self-adjoint operator from K{t) 
to A(t)* □ 

The curve A(-) in L{Ti) is called regular, monotone, monotone increasing (decreasing), 
if its velocity A(t) at any point t is respectively a non-degenerated, sign-definite, positive 
(negative) definite quadratic form on the space A(t). 

Proposition 1 A dynamical distribution {7i,T>) is regular, monotone, monotone increas- 
ing (decreasing) if and only if all Jacobi curves w.r.t. this distribution are respectively reg- 
ular, monotone, monotone increasing (decreasing) curves in the corresponding Lagrange 
Grassmannians. 

Proof. Recall that for any two vector fields H. and £ in M one has 

j^{{e''%e)={e-'%[n,i]. (2.3) 

Let Q^^'^^ be as in (|1.6|) . Applying this fact to the Jacobi curve J\{t) and using (|1.6|) . 
(|2.1() . and (|2.2() one obtains easily that 

= j;,(0), (2.4) 
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which imphes the statement of the proposition. □ 

Fix some A € L{T,). Define the hnear mapping B\ : S i— > A* in the fohowing way: for 
given w € T, one has 

Ba{w){v) = a{w,v), VwgA. (2.5) 

Denote by A^ the set of all Lagrangian subspaces of S transversal to A, i.e. A* = {F € 
Gm{W) : F n A = 0}. Then for any subspace F G A'^ the restriction B\\^ : F i— > A* is an 
isomorphism. 

Remark 5 In other words, any F € A*^ can be canonically identified with the dual space 
A*. 

Let Ir = (-Ba|p) ^ ■ Note that by construction Ir is linear mapping from A* to F and 

Irl - IaI G a. (2.6) 

The crucial observation is that the set A^ can be considered as an affine space over the 
linear space Quad(A*) of all quadratic forms on the space A*. Indeed, one can define the 
operation of subtraction on A* with values in Quad(A*) in the following way: 

{T-A){l) = a{IrljAl). (2.7) 

It is not difficult to show that A*^ endowed with this operation of subtraction satisfies the 
axioms of affine space. For example, let us prove that 

(F - A) + (A - n) = (F - n) (2.8) 

Indeed, using skew-symmetry of a and relation (|2.6j) . one has the following series of iden- 
tities for any / G A* 

(T{Irl, IaI) + (t{IaI, lul) = cj{ItI - IaI) = crihl - IaI - IuL)+ 

a{hl - Inl, Inl) = <y{Irl - lul, M) = <y{hl, IrI), 
which implies (|2.8|1 .^ 

Consider now some curve A(-) in L(S). Fix some parameter r. Assume that A(t) € 
A(r)* for all t from a punctured neighborhood of r. Then we obtain the curve t i— > 
A(t) G A(r)* in the affine space A(r)*. Denote by A^(t) the identical embedding of A{t) 
in the affine space A(r)^. Fixing an "origin" A in A(r) we obtain a vector function 
t I— > Arit) — A with values in Quad (A*). The curve A(-) is called ample at the point r if 
the function Arit) — ^ has the pole at t = r (obviously, this definition does not depend 
on the choice of the "origin" A in A(r)*). In particular, if A(-) is a regular curve in L(S), 
then one can show without difficulties that the function t i— > AT-(t) — A has a simple pole 
at t = T for any A G A(t)'^. Therefore any regular curve in L{T,) is ample at any point. 

^For slightly different description of the affine structure on see and also 4 , where a similar 

construction is given for the Grassmannian Gn(R^") of half-dimensional subspaces of R^". 
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Suppose that the curve A(-) is ample at some point r. Using only the axioms of affine 
space, one can prove easily that there exist a unique subspace A°(r) € such that the 
free term in the expansion of the function t i— > Ar{t) — A°(t) to the Laurent series at r is 
equal to zero. The subspace A°(t) is called the derivative subspace of the curve A(-) at 
the point r. If the curve A(-) is ample at any point, one can consider the curve r ^ A°(r) 
of the derivative subspaces. This curve is called derivative curve of the curve A(-). 

Now assume that the derivative curve A°(t) is smooth at a point r. In particular, the 
derivative curve of a regular curve in the Lagrange Grassmannian is smooth at any point 
(see, for example, |2j and the coordinate representation below). In general, the derivative 
curve of the ample curve is smooth at generic points (the points of so-called constant 
weight, see jSj). As was mentioned already in Remarkl^ one can look on A(r) and A°(t) 
as on the corresponding self-adjoint linear mappings: 

A(r) : A(r) ^ A(r)*, A°(r) : A°(r) ^ (A°(t)) * (2.9) 

Besides, by construction A°(r) G A(r)^. Therefore by RemarklSJthe following spaces can 
be canonically identified: 

A(T)*^A°(r), (A°(t))* ^ A(r) (2.10) 

After these identifications, the composition A°(r) o A(r) is well-defined linear operator on 
A(t). 

Definition 1 The linear operator 

Ra{t) = -A°{t)oA{t). (2.11) 

on A(r) is called the curvature operator of the curve A(-) at a point r. The quadratic form 
^a(''") on A(r), defined by 

rA(r)(i;) = (A(r) o Ra{t)v) (v), v G A(r) (2.12) 

is called the curvature form of the curve A(-) at a point r. 

Suppose that for all Jacobi curves w.r.t. the dynamical Lagrangian distribution (7Y, T>) 

the curvature operator is defined. The curvature operator R^^'^^ of the dynamical La- 
grangian distribution (TC, T>) at a point A is by definition the curvature operator of the 
Jacobi curve J\{t) at t = 0, namely, 

<'^)=i?j,(0). (2.13) 

By construction, it is the linear operator on T>x. Tlie curvature form r^^'^^ of the dynam- 
ical Lagrangian distribution (Ti, T>) at a point A is by definition the curvature form of the 
Jacobi curve Jx{t) at t = (see also (jl.Tf) ). 
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Now for a regular curve A(-) in the Lagrange Grassmannian L{T,) let us give a coordi- 
nate representation of the derivative curve and the curvature operator. One can choose a 
basis in S such that 

S ^ M" X M" = {(x,y) : x,y € M"}, (2.14) 
cr{{xi,yi), (3:2,2/2)) = {xi,y2) - {x2,yi), (2.15) 

where (•, •) is the standard inner product in M" (such basis is called symplectic or Darboux 
basis). Denote by the ith vector of the standard basis of M". 

Assume also that A(t) n {(0,y) : y € M"} = 0. Then for any t sufficiently closed to 
r there exits the symmetric nx n matrix St such that A(t) = {{x,Stx) : x G M"}. The 
matrix curve t 1-^ St is the coordinate representation of the curve A(-) (w.r.t. the chosen 
symplectic basis in Sj. 

Remark 6 Note that from Ii2.1,5]) and the fact that the subspaces A{t) are Lagrangian it 
follows that the matrices St are symmetric. 

The curve A(-) is regular if and only if the matrices St are non-degenerated. The 
expression of the derivative curve and the curvature operator of the regular curve A(-) in 
terms of St is given by the following 

Proposition 2 The derivative curve A°(r) of the regular curve A(r) in L{Ti) satisfies 

A°(^) = {i-\s-^SrS-'y,y- ^SrS;'SrS;'y),y e M"}. (2.16) 

In the basis {{ci, Srei)}'^^^ 0/ A(r) the curvature operator R/\_{t) is represented by the 
following matrix 

S{St) = \s;'si'^-l{S-'Sr)'. (2.17) 
The curvature form rj^{T) has the following matrix w.r.t. the same basis 

~Sr^{St) = -^5(3) + -^SrS;^Sr. (2.18) 

For the proof of (|2.16j) and of the matrix representation (|2.17j) for the curvature operator 
see, for example, jl]. The matrix representation 1)2. 18(1 of the curvature form follows 
directly from (|2T7|) and ((2?T2|) . 

Remark 7 If St is a scalar function (i.e., n = 1), then S{St) is just the classical Schwarzian 
derivative or Schwarzian of St. It is well known that for scalar functions the Schwarzian 
satisfies the following remarkable identity: 



§{ip{t)) (2.19) 



\cip{t) + 

for any constant a, b, c, and d, ad — be ^ 0. Note that by choosing another symplectic 
basis in S, we obtain a new coordinate representation i i— > jSj of the curve A(-) which is a 
matrix Mobious transformation of St, 

St = iC + DSt)iA + BSt)-^ (2.20) 
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for some n x n matrix A, B, C, and D. It turns out that the matrix Schwarzian ()2.17j) 
is invariant w.r.t. matrix Mobious transformations (|2.'20|) by analogy with identity ()2.19|) 
(the only difference is that instead of identity we obtain similarity of corresponding ma- 
trices). This is another explanation for invariant meaning of the expression H2.17() . given 
by Proposition |2I 

The coordinate representations (|2.17j) and (|2.18|) are crucial in the proof of the main 
theorem of this section (see Theorem ^ below). 

2.2 Curvature operator and curvature form of reduction. Now fix some s 
vectors /i, . . . in S such that 

yi<i,j<s: a{li,lj) = 0. (2.21) 
Denote by ^ = {h, ■ ■ ■ h) and span^ = span(/i, . . .Is) For any A G L(Y,) let 

A^ = An (span £)^ + span F = AVspan^, (2.22) 

where (span£)^ {w € S : VI < z < s a{v, li) = 0} is the skew-orthogonal complement 
of the isotropic subspace span i. Actually, A^ is a Lagrangian subspace of the symplectic 
space (span^)^/span£ (with symplectic form induced by a). 

Let, as before, A(-) be an ample curve in L{Ti). The curve A(-)^ is called the reduction 
by the s-tuple £, satisfying 1)2. 21() . or shortly the i-reduction of the curve A(-). Note 
that by (|2.22|) span£ C A(f)^ for any t. Therefore the curve A(-)^ is not ample and the 
constructions of the previous subsection cannot be applied to it directly. Instead, suppose 
that the curve A(-)^ is ample curve in the Lagrange Grassmannian L((span£)^/span£). 
Then the curvature operator R-j^i^t) of this curve is well-defined linear operator on the 

space A{tY (at least for a generic point t). Let cp : T, S/span^ be the canonical 
projection on the factor-space. 

Definition 2 The curvature operator R^e{T) of the (-reduction A(-)^ at a point r is the 
linear operator on A(tY , satisfying 

RMir)iv) = ('AL(.)n(span.)-)"' ° ^af(^) ° '^(^)' ^ e A{tY. (2.23) 

The curvature form r^e{T) of the i-reduction A(-)^ at a point r is the quadratic form on 
A(t)^, satisfying 

rAr){v) = ^{K{tY){R^,{t)v,v). (2.24) 

All these constructions are directly related to the reduction of dynamical distributions 
by a tuple Q = (gi, . . . ,gs) of s involutive first integrals, defined in Introduction. Indeed, 
the Jacobi curves attached at some point A w.r.t. the ^-reduction {Tl,T>^) of a dynami- 
cal Lagrangian distribution {H,T)) are exactly (^i(A), . . . , 5s( A)) -reductions of the Jacobi 

curves attached at A w.r.t. (7i, V) itself. The curvature operator R^^'^ ^ and the curva- 
ture form r^x^''^ ^ at X of the Q-reduction (7Y, D^) are by definition the curvature operator 
and the curvature form of the Jacobi curves attached at A w.r.t. {TC,T>^). 
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To justify these definitions suppose that we are in situation of ExamplelHl i.e. 7i admits 
one first integral g, satisfying (|1.15|) . Let a symplectic manifold Wg^c and a mapping 
$ : g~^{c) ^ Wg^c be as in this example. Then directly from the definition it follows that 

V^eP^: rf^'''\v)=rf^f''*'''\<^>.v). (2.25) 

In other words, the curvature form of the ^(-reduction {H,V^) at A G 9~^{c) is equal 
to the pull-back by ^> of the curvature form of the dynamical Lagrangian distribution 
{^jfTtj^^T)^), associated to the original dynamical Lagrangian distribution {h,T>) on the 
reduced symplectic space Wg^c- 

The natural question is what is the relation between the curvature forms and operators 
of the dynamical Lagrangian distribution and its reduction on the common space of their 
definition. Before answering this question in the general situation, let us consider the 
following simple example: 

Example 8 (Kepler's problem) Consider a natural mechanical system on M = with 
the potential energy U = —r~^, where r is the distance between a moving point in a plane 
and some fixed point. This system describes the motion of the center of masses of two 
gravitationally interacting bodies in the plane of their motion (see ) . Let q = (r, if) 
be the polar coordinates in M^. Then the Hamiltonian function of the problem takes the 
form 

2 2r2 r 

where pr and are the canonical impulses conjugated to r and ip. If A = {p,q), where 
q e M, p £ T*M, then pr{\) = p{dr{q)) = dgV, p^{X) = r'^p{d^{q)) = r'^dqtp. Observe 
that g = p^ is nothing but the angular momentum and from (|2.26|) we immediately see 

that it is a first integral of the system. Let us compare the curvature forms r^^'^^^'^')^ and 
^(/i.,n(M) ) ^j^^ common space n(M) n ker dxg = Mc?p^ of their definition. 

First, according to (|1.9() of Example 2, the curvature form of (/i, n(M)) is equal the 
Hessian of U at q. In particular, it implies that 

(h,n(M)) , 9^ , , 2 , , 

rl '\d,^) = —^U{q) = --^. (2.27) 

Further, let c = g{X). Note that g satisfies the condition H1.16() of Remark|21with V = 
r'^d^. Let Wg^c and $ be as in Example El Then, following Remark 01 Wg^c = r*M+ and 
the dynamical Lagrangian distribution (<I>*/i, <I>^,M9p^) is equivalent (symplectomorphic) 
to the dynamical Lagrangian distribution associated with the natural mechanical system 
with configuration space M"*" and the potential energy 

u 

2r^ r 

(Ua is the so - called amended potential energy; it comes from the following identity: 
^L-i(c) = f + Uair)). Hence by 



^ = ?f + S-r> (2-26) 



ri'^'-''^ iOp.) - ^U.ir) ^ ^ - ^ ^ .(-^-)) (..J + ^. (2.28) 
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Note that from H2.28() it follows that on the common space of the definition the reduced 
curvature form is not less than the curvature form itself. We will show later (Corollary |2) 
that this is a general fact.D 

2.3 The change of the curvature after the reduction. Now we give the relation 
between the curvature forms of the regular curve A(-) and its ^-reduction A(-)^, where ,as 
before, i = {h,. . . ,ls) is the tuple of s vectors, satisfying (|2.2H) . First, let us introduce 
some notations. Let -BA(t) : S i— > A* be as in (|2.5() . Looking at A{t) as at a linear mapping 
from A(t) to A(t)*, denote by aj(t), 1 < i < s the following vectors in A(t): 

aiit) = {A{t))-'oB^^t){h). (2.29) 

Using definition of A{t) and i?A(t) can show that t i— ai{t), aiit) € A(t), is a unique 
vector function such that for any t one has 

a{ai{t),v) = a{li, v), Vw G A(t), 

or, equivalently, 

ai{t) = k mod A{t). (2.30) 
Finally, let A{t) be the s x s matrix with the following entries 

Akm{t) = cF{lk,am{t)), l<k,m<s. (2-31) 

Note that by (|On|) and the definitions of A{t) 

Akmir) = (j{ak{t),a^{t)) = A{t)ak{t){a^{t)) , (2.32) 

which implies that the matrix A{t) is symmetric. 

Theorem 1 Suppose that A(-) is a regular curve in the Lagrange Grassmannian L{Ti) and 
i = (li, . . . ,ls) is a tuple of s vectors in T, such that \2.21]) holds and det A{t) ^ for some 
point T. Then the curvature form r\{T) of the curve A(-) and the curvature form r/^i{t) of 
its £ -reduction A{-y at the point t satisfy the following identity for all v € A(r)n (span£)^; 



rj,e{r){v) - rA{r){v) = ^ E {Mr)-') {dkir) , v) a {dm{r) , v) , (2.33) 

k,m=l 

where (A(r)^^)^^ is the km-entry of the matrix A{t)~^ . In addition, the curvature oper- 
ator Ra{t) of the curve A(-) and the curvature operator R^e{T) of its i-reduction A(-)' at 
the point r satisfy on A(r) R (span£) the following identity: 

3 

i?A.(r) - Ra{t) = ^Y1 iMrr')kmBA^r)amir)(^[{HT)y' o B^^^dkir)) . (2.34) 

k,m=l 

(As usual, for a given linear functional ^ and a given vector v by ^ <Si v we denote the 
following rank 1 linear operator v{-) = £,{-)v.) 
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Proof. First let us prove identity (|2.33() . As before, denote by Cj the ith. vector of the 
standard basis of M". The condition det j4(t) ^ it equivalent to 

span(ai(r), . . . , asir)) n (span^)^ = 0. (2.35) 

Hence one can choose a basis in S such that if one coordinatizes S w.r.t. this basis, 
E = X M", then the symplectic form a satisfies (|2.15|) and the following relations hold 

A(r) n (span^) ^ = span ((ei, 0), . . . , (e^-., 0)) , (2.36) 
li = {0,en-s+i), l<i<s, (2.37) 
ai(T) e span((en-s+i,0), . . . , (en,0)), l<i<s. (2.38) 

Note that by construction 

(span^)^ = span ((ei,0), . . . , (e„_s,0), (0, ei), , . . . , (0, e„)). 

Therefore one can make the following identification: 

(span^)^/span£ = span ((ei,0), . . . , (e„_s,0), (0,ei), . . . , (0, e„_s)). (2.39) 

Since by definition ai{t) e A(t), there exists bi{t) € M" such that a{t) = (bi{t), Stbi{t)) . 
Note that 

ai{t) = {0,Stbi{t)) + {bi{t),Stbi{t)) = {0,Stbi{t)) modA(f). (2.40) 
This, together with (|On|) and ((233), implies that k = {0,Stbi{t)) and then 

bi{t) = S^'cn-s+i. (2.41) 
On the other hand, from (|2.37|) and (|2.38j) . using i|2.15j) one can obtain that 

s s 

bi{T) = '^a{ai{T),lj)en-s+j = - ^ ^ij(r)e„_s+j (2.42) 
i=i j=i 

(in the last equality we used the symmetry of the matrix A{t)). So, from 1)2. 41(1 . (|2.42|) 
and symmetry of Sr (see Remark El it follows that 

yi <i<n-s, n-s + l<j<n: (5'^i),y = (57^)^^ = 0, (2.43) 
yi<i,j<s: {S;%^s+i,n-s+j = -Aijir). (2.44) 

Further, for given nxn matrix A denote by C{A) the {n — s)x (n — s) matrix, obtained 
from A by erasing the last s columns and rows. Consider the curve A(-)^ in the Lagrange 
Grassmannian L((span£)^/span£) (see (|2.22|) for the notation). By construction, if St is 
the coordinate representation of the curve A(t) w.r.t. the chosen symplectic basis, then 
C{St) is a coordinate representation of the curve A(-)' w.r.t. the basis of (span£)^/span^, 
indicated in (|2.39jl . Hence from (|2.43jl . (|2.44j> . and assumption det A{t) 7^ it follows 



16 



that the germ at r of the curve A(-)^ is regular. In particular, the curvature form of the 
^-reduction A(-)^ is well defined at r. Using H2.18() . we obtain that the quadratic form 

A(r)n(spanf)^ A(T)n(span 

has the following matrix in the basis ((ei, 0), . . . , (e^-s, 0)) : 

-^C{Sr)S{C{Sr)) +C{SrS{Sr)). 

Using the blocked structure of the matrix Sr, given by H2.43() . one can obtain from ()2.18|) 
without difficulties that 

/7 -if Tl — S 

- — C{Sr)S{CiSr)) + C{SrS{Sr)) = --{ {SrMS-\miSr)mj } . .^^ 

k,m=n-~s+l 

In order to prove (|2..S3j) . it is sufficient to prove the following 
Lemma 1 The restriction of the quadratic form 



s 

V I— > 

k,m=l 



(.A{t) ^)kni(T{ak{T),v)a{am(,T),v) (2.45) 

t,m=l 

on A(r) n has the matrix with ij -entry equal to 

n 

~ ^ ^ {S-r)ik{Sj- )km{S-r)mj 
k,m=n—s+l 

in the basis ((ei, 0), . . . , (e^-s, 0)) . 

Proof. First, using the symmetry of St and (|2.44|) . one has the following identity: 

n n 
^ ^ ('S'T)ifc('S'T- )km{ST)mj = ^ ^ {SrSj. )j,fc(5'T)fcm('5r5'T- )jm = 

• ..''"7""' , .. . , 

k,m=l 

On the other hand, from (|2.4()|) we have: 

■diit) = {0,Stbi{t)) +2{0, Stk{t)) + (bi{t),SMt)). 

Substitute i = r in the last relation. Note that Sr has the same blocked structure, as 
S~^, which together with (|2.42|) imphes that (0,Stbi{t)) € span£. Prom this and (|2.41|1 it 
follows that 

diir) = {0,SrS~^en-s+iit)) mod span(A(r), /i, . . . , Z^) . (2.47) 
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Hence 

\Jl<j<n-s,l<i<s: a((di{T),{ej,^))=-{SrS;^)j,n-s+u (2.48) 
The last identity together with (|2.4(i)l imphes the statement of the lemma and also formula 

Finally, identity (|2.l-i4|) follows directly from (|2.12|) . The proof of the theorem is com- 
pleted. □ 

Note that if the curve A(.) is monotone increasing or decreasing , then from (|2.32|) it 
follows that the condition det A{t) 7^ is equivalent to the following condition 

A(t) nspan^ = 0. (2.49) 

As a direct consequence of identity (|2.33() and the last fact one has the following 

Corollary 1 // the curve A(.) is monotone increasing and the tuple £ = {li, . . . ,ls) of s 
vectors in S satisfies \2.4y^ for some t, then on the space A{t) n (span^)^ the curvature 
form of the i-reduction of the curve A{-) is not less than the curvature form of the curve 
A(.) itself. Moreover, on the space A(t) n (span^) the difference between the curvature 
form of the i-reduction of the curve A(.) and the curvature form of the curve A(.) itself is 
non-negative definite quadratic form of rank not greater than s. 

Let us translate the results of Theorem ^ in terms of a regular dynamical Lagrangian 
distribution {H, V) on a symplectic space W . Suppose that the Hamiltonian Ti admits 
a s-tuple Q = (gi, . . . ,gs) of involutive first integrals. Similarly to H2.5() denote by Bd^ '■ 
T\W ^ 2?^ the linear mapping such that for given Y ^T\W the following identity holds 

Bv^Y{Z) = a{Y, Z), VZ G V^. (2.50) 

Let us look at Ja(0) (the velocity at t = of the Jacobi curve attached at A) as at a linear 
mapping from "Dx to P^. Then using regularity by analogy with 1)2. 29|) one can define the 
following s vector fields Xi on W: 

Xi{\) = {jx{0)y^oBv,{gi{X)). (2.51) 

Using relation (|2.3j) one can obtain by analogy with (|2.3()|) that Xi is a unique vector field, 
satisfying Xi{X) € Vx and 

[n,X^]{X)=g,{\) mod Pa (2.52) 
for all A G W. Finally let T(A) be the s x s matrix with the following entries 

'^Wkm = crx{gk,Xm), l<k,m<s. (2.53) 
Prom Theorem Hand relation (|2.3|) one gets immediately 
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Theorem 2 Suppose that {Ti.,T>) is a regular Lagrangian dynamical distribution on a 
symplectic space W and the Hamiltonian TL admits a tuple Q = {gi, . . . , gs) of s involutive 

first integrals such that det T(A) 7^ 0. Then the curvature form r^^'^^ of the dynamical 

distribution (??, P) and the curvature form r^^'^ ^ of its Q-reduction {7i,T>^) satisfy the 

s 

following identity for all v ^ \ {^ d\gi j H T>\ 



i=l 



fc,m=l 

(2.54) 

while the curvature operator R^^'^^ of the dynamical distribution (Ti, V) and the curvature 

s 

operator R^^'^ ■* of its Q-reduction {'H,'D^) satisfy on ker d\gij fl Vx the following 

i=l 

identity: 

^ ^ (T(A)-i),^^?»j7?,[7i,^^]](A)®((j,(0))-^S^j7?,[7?,^fc]](A)). ^^'^^^ 



4 

k,m=l 

Also, by analogy with Corollary^ we have 

Corollary 2 // the dynamical Lagrangian distribution (7Y, V) is monotone increasing and 
the Hamiltonian Ti admits a tuple G = {gi, ■ ■ ■ , gs) of s involutive first integrals such that 

DAnspan(5i(A),...,5s(A)) =0, (2.56) 



then on the space ker d\gij the curvature form of the Q-reduction of the dynam- 

i=l 

ical Lagrangian distribution (Ti.,!)) is not less than the curvature form of (?{,T>) itself. 

s 

Moreover, on the space ^f^ker dxgi^ ri'Dx the difference 

1=1 

'^A '^A 

is non-negative definite quadratic form of rank not greater than s. 

Now let us give the coordinate representation of the vector fields Xi, \ < i < s from 
Theorem El in the case, when W = T*M and V = n(M). Let q = (g\ . . . ,g") be local 
coordinates in some open subset M oi M and p = (pi, . . . ,p„) be induced coordinates 

n 

in the fiber of T*J\[ so that the canonical symplectic form is given by cj = ^ dpi A dq^. 

i=l 
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It gives the identification of T*M = {{p,q),p,q G M") (so, J\f = xW). 

Also the tangent space TxiT*M) to T*M at any A is identified with M" x M". Under 

this identification H = (^~^)^^) where for given function h on T*Af we denote by 

i = (1^' • • • ' and f = (1^, . . . , 1^)^. Denote by H^p the Hessian matrix of 

the restriction of Ti to the fibers. Then from Remark ^ and relation (|2.5H) we have 

^^ = iKp^^O) (2-57) 

Now suppose for simphcity that the dynamical Lagrangian distribution is associated 
with a natural mechanical system (Example 2) or, more generally, with a mechanical 
system on a Riemannian manifolds (Example 3) , which admits one or several first integrals 
being in involution and linear w.r.t. the impulses. One way to compute the reduced 
curvatures is to pass to the reduced phase space, as was described in Remark |31 and apply 
the method of computation of the curvatures from |^ to the corresponding dynamical 
Lagrangian distribution in the reduced phase space (this way was actually implemented 
in Example ISJ. But in order to apply the method of flj we need to find a new canonical 
coordinates in the reduced phase space, which is not just a trivial exercise. Moreover, very 
often the new Hamiltonian system on the reduced phase space has more complicated form 
than the original one. Both these facts make the computation in this way quite tricky. 
Theorem 2 gives another method to compute all reduced curvatures without passing to the 
reduced phase space: to do this one can combine (|1.9|) or (|1.12|) with (|2.54|) (or (|2.55j) ) and 
1)2. 57(1 . This method is more effective from the computational point of view, especially if 
the number of the involutive first integrals is essentially less than the number of the degrees 
of freedom in the problem. We illustrate the effectiveness of this method on the following 
example: 

Example 9 ( Plane N-body problem with equal masses) Let us consider the motion of 
bodies of unit mass in endowed with the standard Cartesian coordinates so that 
= {q2i-i,Q2i) G represents the radius vector of the i-th. body with respect to some 

inertial frame. It is described by a natural mechanical system on M = with potential 

energy 

^ 1 

C/(ri, . . . ,rAr) = — , nj = \\ri - rjW . (2.58) 
Then T*M ^ M^TV ^ = {(^^ g)^ p^q(z ^2^}^ 

Pi ■ ■ ■ , P2N are the canonical impulses 
conjugated to gi, . . . , q2N (Pi ^ Qi)- The systems has the following first integral 

N 

g = '^{P2iq2i-i - P2i-iq2i) (2.59) 
1=1 

which is nothing but the angular momentum (in the considered planar case the angular 
momentum is scalar). From Example 2 we know that the generalized curvature form of the 
dynamical Lagrangian distribution (i?,n(M)) is just the Hessian of the potential energy 
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U and the generalized Ricci curvature (see ()1.8j) for the definition) is the Laplacian of U, 
which can be calculated without difficulties: 



Px 



{h,U{M)) 



^ 1 



i<j v 



(2.60) 



Our goal is to compute the reduced generalized Ricci curvatures ' ^ , using 
the formula (|2.54|) . In our case s = 1. Let X be as in (|2.5H) with g instead of gi. Note 
that by definition the vector -^(A) is orthogonal to the subspace Il{M)x n kerd^^^f w.r.t. 

the mner product I 



•, •). Therefore 



Px 



tr 



H,n{M)3) (/?,n(M)) 



(j?,n(M)) 



R 



Px 



(2.61) 



A 



R 



'h,II{M)) 



n(Af)AnkcrdA9. 



Further, the Hamiltonian vector fields corresponding to the functions H and g are given 

hy ii = {Uq,p) and g = {Jp,—Jq) with Uq = and J being the unit 

symplectic 2N x 2N matrix: 

/ 1 ... \ 
-10 ... 

J = 

... 1 

V ... -10/ 

Applying formula (|2.57|1 we find that 



X = iJq,0). 



(2.62) 



Denote X = Jq. Also, let (■,•) and || • || be the standard Euclidean inner product and 
norm. Using again ()1.9|) . one can obtain by direct computation that 



'^x 



iX) = {UqqX,X) 



-u. 



(2.63) 



Further, using (|2.54|) (or (|2.55j) ') and (|2.62|) . one can obtain without difficulties that 



tr 



R 



(i?,n(M)9) (h,u(m)) 

— Ux 



Aa{X,g) 



2N 



STrrdn iff y]] ri o-([-g, [H,X]],xy 



i=l 



(2.64) 
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Substituting dTHHl) . fT^ . and dTHl) into itTHTl) we obtain finally that 



.P*'"'=-2E;i-T + ^P"4wn^ (2.65) 

where / =||g|p,T=^||p|p are the central momentum of inertia and the kinetic energy of 
the system of N bodies. Note that the sum of the first two terms in (|2.65|) is the trace of 

the restriction of the curvature operator ^ on the space Il{M)\ n kerdA^. So, by 

Remark[2land CorollarylJl the last term in (|2.65|) has to be nonnegative. Actually this term 
contains the right-hand side of the famous Sundman's inequality 2TI — \{H, I}^ > and 
it is nothing but the generalized area of the parallelogram formed by two 2A^-dimensional 
vectors p and q. □ 



3 Focal points and Reduction 



In the present section we study the relation between the set of focal points to the given 
point w.r.t. the monotone increasing (or decreasing) dynamical Lagrangian distribution 
and the set of focal points w.r.t. its reduction. As before, first we prove the corresponding 
result for the curves in Lagrange Grassmannians and then reformulate it in terns of the 
dynamical Lagrangian distributions. 

Let A(-) be a curve in the Lagrange Grassmannian defined on the interval [0, T]. 

The time ti is called focal to the time w.r.t. the curve A(-), if A(ti) n A(0) 7^ 0. The 
dimension of the space A(ti) n A(0) is called the multiphcity of the focal time ti. Denote 
by 7^focoA(-)|^ the number of focal times to on the subset I w.r.t. A(-), counted with 
their multiplicity. If the curve A(-) is monotone increasing, then 7^focoA(-)|^ is finite and 
one can write 

#focoA(-)|, = J^dim(A(t) n A(0)). (3.1) 

tei 

Fix some tuple i = {h. . . ,ls) of s linearly independent vectors in E, satisfying 1)2.21(1 . 
The time ti is called focal to the time w.r.t. the £-reduction A(-)^ of the curve A(-), if 
A^(ii) n A^(0) 7^ span£ or, equivalently, ti is the focal time to w.r.t. the curve A(-)^ in 
the Lagrange Grassmannian L((span£)^/span£). The multiplicity of the focal time ti to 
w.r.t. the ^-reduction A(-)^ is equal by definition to 

dim {MhY n A(0)^) = dim {A{tiY n A(0)0 - s. (3.2) 

So, the number of focal times to on the subset / w.r.t. the .^-reduction A(-)^, counted 
with their multiplicity, is equal to #focoA(-)^|^. 

It is not hard to see that if A(-) is monotone increasing, then A(-)^ is monotone in- 
creasing too. So, the number of focal times to w.r.t. the ^-reduction is also finite. The 
question is what is the relation between the set of points, which are focal to w.r.t. the 
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curves A(-) and its ^-reduction A(-)^? To answer this question, we use the fact that if A(-) 
is monotone curve, then the number #focA(-)|^Q can be represented as the intersection 
index of this curve with a certain cooriented hypersurface in L{T,). The advantage of this 
representation is that the intersection index is a homotopic invariant. 

More precisely, for given Lagrangian subspace Aq denote by A^Aq the following subset 

of 

Ma, = L(S)\A* = {A G : A n Aq / 0} . 

Following [HI, the set A^Aq is called the train of the Lagrangian subspace Aq. The set 
A^Ao is a hypersurface in L{T,) with singularities, consisting of the Lagrangian subspaces 
A such that dim (An Aq) > 2. The set of singular points has codimension 3 in L{T,). As we 
have already seen, the tangent space Ta-L(S) has a natural identification with the space 
of quadratic forms on A. If A is a non-singular point of the train Maq , then vectors from 
TaL(Y,) that correspond to positive or negative definite quadratic forms are not tangent to 
the train. It defines the canonical coorientation of the hyper-surface A4ao a-t a non-singular 
point A by taking as a positive side the side of Maq containing positive definite forms. 
The defined coorientation permits to define correctly the intersection index A(-) ■ A^Aq 
of an arbitrary continuous curve in the Lagrangian Grassmannian A(-), having endpoints 
outside A^Ao- If ^(O is smooth and transversally intersecting A^Aq i^i non-singular points, 
then, as usual, every intersection point A{t) with A^Aq adds or —1 into the value 
of the intersection index according to the direction of the vector A{t) respectively to 
the positive or negative side of A^Aq- Further, an arbitrary continuous curve A(-) with 
endpoints outside A^Ao can be (homotopically) perturbed to a curve which is smooth and 
transversally intersects A^Aq i^i non-singular points. Since the set of singular points of 
A^Ao has codimension 3 in L{T,), any two curves, obtained by such perturbation, can 
be deformed one to another by homotopy, which avoids the singularities of AIaq • Hence 
the intersection index of the curve, obtained by the perturbation, does not depend on 
the perturbation and can be taken as the intersection index of the original curve. The 
intersection index A(-) • A^Aq can be calculated using the notion of the Maslov index of 
the triple of Lagrangian subspaces (see [S] , [El ) for details) . This implies in particular that 
if the curve A : [0, T] i-^ L{T,) is monotone increasing with endpoints outside A^Aq! then 

M-)-Mao= dim(A(t)nAo) (3.3) 

0<t<T 

Now we are ready to formulate the main result of this section: 

Theorem 3 Let A : [0,T] ^ L(Y,) be a monotone increasing curve and I = {li . . . ,ls) be 
a tuple of s linearly independent vectors in T,, satisfying h2.21\) and {2.4-9^ at 0. Then on 
the set (0, T] the difference between the number of focal times to w.r.t. the i-reduction 
A{-Y and the number of focal times to w.r.t. the curve A(-) itself , counted with their 
multiplicity, is nonnegative and does not exceed s, namely 



< #focoA(-)^|(o_^] - #focoA(-)|(o^y] < s . (3.4) 
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Proof. Since the curves A(-) and A(-)^ are monotone increasing, for sufficiently small 
e > the set (0, e] does not contain the times focal to w.r.t. both of these curves. Also, 
without loss of generality, one can assume that T is not focal to w.r.t. both of these 
curves (otherwise, one can extend A(-) as a monotone increasing curve to a slightly bigger 
interval [0,r + e] such that A(r + e) n A(0) = and A(r + eY n A(0)^ = 0). So, by 
relations (EH), (ESI), we have first that 

#focoA(-)|(o,r] = ^(Olte.T] ■ ^m- (3-5) 
Besides, from (|2.49j) it follows that 

dim {A{tY n A(0)^) - s = dim (A^(t) n A(0)) . (3.6) 
Hence, combining (|3.H) and (|.S.3|) with (|3.2|) and (|3.6|) . we get 

#focoA(7|(Q^^j = A(.)^|[^_^] • >tA(o) (3.7) 

Now we prove the theorem in the case s = 1. In this case i = li. We use the invariance 
of the defined intersection index under homotopies, preserving the endpoints. 
Let ai{t) be as in (|2.29|) . Denote 

F(r, t) = span(A(t) n if, (1 - r)ai(t) + r/i) (3.8) 

Note that all subspaces F{t, t) are Lagrangian. Let '■ [0, T] LCS) and Tt : [0, 1] i— >• 
L(S) be the curves, satisfying 

1>.(-)=F(r,-), 0<T<1; 
Tti) = Fi;t), 0<t<T. 

Then $o(-) = A(-), $i(-) = A'i(-), and the curves 

r.(-)|[o,.]U<i>.(-)|[,,^]u(-rT(-)|[o,.]) 

define the homotopy between A(-)|j^^j and re(-) U -^(O'^lieT] ^ {~^t{-)), preserving the 
endpoints (here — 7(-) means the curve, obtained from a curve 7(-) by inverting the orien- 
tation). Therefore, 

^(■)|[e,r] • ^MO) = r.(-) ■ -^A(O) + M-y'\[e,T] ■ -^A(O) " ^t{-) ' >Ja(0) 
Using (|3.5|) and (|3.7j) . the last relation can be rewritten in the following form 

#iocoAiY'\^o,T] - #focoA(-)|(o,T] = rT(-) • Ma^o) - r.(-) • 7Wa(o). (3.10) 

So, in order to prove the theorem in the considered case it is sufficient to prove the following 
two relations: 

< rT(-) • A^A(O) < 1, (3.11) 
3eo > s. t. Veo > e > : T^i-) ■ A4a(o) = 0. (3.12) 
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a) Let us prove XTW\ . If h G A(r), then by definition Tt{t) = A{T). Since, by our 
assumptions, A(T) n A(0) = 0, we obviously have Tt{-) ■ ^^^(o) = 0- 

If h A(r), then dim (A(0) + A(r) n (Zi)^) = 2n - 1. In particular, it implies that 

< dim {TTir) n A(0)) < 1. (3.13) 

Further, let p : S i— > S/(A(0) + A(T) n (h)^) be the canonical projection on the factor 
space. Then from (|3.8() and 1)3. 9() . using standard arguments of Linear Algebra, it follows 
that rr(r) n A(0) 7^ if and only if 

(1 - T)p{aiiT)) + Tp{h) = (3.14) 

Since, by assumptions, Tt{0) n A(0) = (recaU that Tt{0) = A{T)), the equation (fXTH) 
has at most one solution on the segment [0, 1]. In other words, the curve T -];•{•) intersects 
the train A1a(o) most ones and according to 1)3. 13(1 the point of intersection is non- 
singular. 

Finally, the curve rr(-) is monotone non-decreasing, i.e. its velocities ^r7^(r) are 
non-negative definite quadratic forms for any r. Indeed, since A(r) n (h)'^ is the common 
space for all Tt{t), one has ^rr(''")|A(T)n(/i)'^ ^ ^' other hand, if we denote by 

c(r) = (1 — T)ai{t) + r/i, then by (|2.2|) . one has 

^rT(r)(c(T)) =a{c'{T),c{T)) =a{h-ai{T),{l-T)aiiT)+Th) =a{h,ai)>0 

(the last inequality follows from ()2.31|) . ()2.32|) and the assumption about monotonicity of 
A(-))- 

So, ^rT(T) are non- negative definite quadratic forms. Hence in the uniquely possible 
point of intersection of Tt{-) with the train A^a(o) the intersection index becomes equal 
to 1. This proves (|3.11|) . 

b) Let us prove (|3.12j) . Take a Lagrangian subspace A such that Zi G A and AnA(O) = 
0. Then there exists eo such that 

Similarly to the arguments in a), for any < e < eo the curve Ti;{-) intersects the train A4a 
once. But by construction this unique intersection occurs at r = 1. Indeed, r£(l) = A(e)^i , 
hence h G re(l) n A. In other words, 

r.(-)|[o,i)CA*. (3.16) 

Further, one can choose a symplectic basis in S such that S = M" x M", the symplectic 
form a is as in (I2T5I1 . A(0) = x R", and A = M" x 0. By (|3T5|l and (ITTfill . there exists 
two one parametric families of symmetric matrices St, < i < eo and Cr ,0 < r < 1 such 
that A(t) = {iStp,p) : p G M"} and T^{t) = {{CrP,p) : p G M""}. Since the curve A(-) is 
monotone increasing and the curve T^{-) is monotone nondecreasing, for any < r < 1 
the quadratic forms p 1— > {CrP,p) are positive definite, while 5*0 = 0. It implies that 

VO < r < 1 : r,(r) fl A(0) = 0. (3.17) 
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Note also that for sufficiently small e > 

r,(i)nA(o) = 0. (3.18) 

Indeed, T^^l) = A'^ (e) and a sufficiently small e > is not a focal time for the /i-reduction 
A'i(-), which according to (|3.6j) is equivalent to the fact that A'i(e)nA(0) = and hence to 
H3.18() . By ()3.17|) and H3.18() . for sufficiently small e > the curve re(-) does not intersect 
the train A1a(o) • The relation (|3.12|) is proved, which completes the proof of our theorem 
in the case s = 1. 

The case of arbitrary s can be obtained immediately from the case s = 1 by induction, 
using the fact that 

Remark 8 Note that in the case s = 1 from Theorem El it follows immediately that the 
sets of focal times ( to 0) w.r.t monotone increasing curve and its reduction are alternating. 
Also, for any s the first focal time to w.r.t. the reduction does exceed the ffist focal time 
w.r.t. the curve itself. □ 

All constructions above are directly related to the notion of focal points of a dynamical 
Lagrangian distributions and {H^T)) and its reduction by a tuple Q = {gi, . . . , Qs) of s 
involutive first integrals, defined in Introduction. Note that the point Ai = e*^^Ao is focal 
to Aq w.r.t. the pair {Tl, V) along the integral curve 1 1— > e*^Ao of TC if and only if the time 
ti is focal to w.r.t the Jacobi curve Jao(') attached at the point Aq, while Ai = e*^^Ao is 
focal to Ao w.r.t. the ^-reduction of the pair {TL, V) along the integral curve t l—^ e*^Ao of 
H if and only if ti is focal to w.r.t {gi{Xo), ■ ■ ■ , ft9(Ao))-reductions of the Jacobi curves 
Jao(-) attached at Aq. Translating Theorem |21 into the terms of dynamical Lagrangian 
distribution, we have immediately the following 

Corollary 3 // the dynamical Lagrangian distribution (Ti, V) is monotone increasing and 
the Hamiltonian TL admits a tuple Q = (gi, . . . , gg) of s involutive first integrals satisfying 
112.56]) . then along any segment of the integral curve ofTL the difference between the number 
of the focal points to the starting point of the segment w.r.t. the Q -reduction of the pair 
(TL,T>) and the number of the focal points to the starting point of the segment w.r.t. the 
pair (TL,T>) itself , counted with their multiplicity ^, is nonnegative and does not exceed s. 

Example 10 (Plane 3-body problem with equal masses: focal points of the 8-shaped 
orbit) The following example illustrates the Theorem |31 In 2000, A. Chenciner and 
R.Montgomery proved the existence of a new periodic solution of the planar 3-body prob- 
lem with equal masses - the 8-shaped orbit or just the Eight |9j. In the plane of the 
motion each body moves along the same 8-shaped orbit, symmetric w.r.t. the point of its 
self-intersection, coinciding with the center of mass of the bodies. The configuration space 
is M = M^. As an initial point in the phase space T*M we take the point Aq such that 

''Here we do not count the starting point of the segment as the focal point to itself. 
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its projection on the configuration space is a colhnear configuration, i.e. one of the bodies 
hes in the middle of the segment, connecting the other two. 

In [H] there were found the focal points to Aq along the Eight w.r.t. the g'-reduction 
of the Lagrangian dynamical distribution (i?, n(M)3), where H, g are as in Example 9. 
In particular it was shown numerically that the 8-shaped orbit contains three such focal 
points along its period T, and the first focal time ri « 0.52T. 

Let e^'^Xo be the ith focal point w.r.t. [H, n(M)) along the Eight, and let e'^'^Xo 
be the ith focal point w.r.t. its g-reduction n(M)^) along the same curve. In the 
following table we present the result of the numerical computation of tj and Tj on the 
interval (0,3T] (all this focal points have the multiplicity 1): 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 




0.52 


0.76 


0.95 


1.08 


1.52 


1.56 


1.88 


2.05 


2.29 


2.49 


2.65 


ti/T^ 


0.76 


0.95 


LOS 


1.42 


1.54 


1.88 


2.05 


2.28 


2.45 


2.65 





We observe that Tj < < Tj+i, 1 < i < 10, as was expected by Theorem 01 □ 
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